First principles investigations of the electronic, magnetic and chemical bonding 

properties of CeTSn (T = Rh, Ru) 
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The electronic structures of CeRhSn and CeRuSn are self-consistently calculated within den- 
sity functional theory using the local spin density approximation for exchange and correlation. In 
agreement with experimental findings, the analyses of the electronic structures and of the chemical 
bonding properties point to the absence of magnetization within the mixed valent Rh based system 
while a finite magnetic moment is observed for trivalent cerium within the Ru-based stannide, which 
contains both trivalent and intermediate valent Ce. 

PACS numbers: 07.55.Jg, 71.20.-b, 71.23 
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I. INTRODUCTION 



Equiatomic cerium-transition metal (T)-stannides 
CeTSn have intensively been studied over the past twenty 
years with respect to their outstanding magnetic prop- 
erties. CeNiSuji CeRhSnfS, and CelrSn'^ are intermedi- 
ate valence systems, while CePdSn (Tjv = 7.5 K)^ and 
CePtSn {Tn = 8 K) (see Ref. H and references therein) 
are antiferromagnetic Kondo lattices. CeAgSn orders an- 
tiferromagnetically at 6.5 and CeAuSn (Tc = 4.1 K)^ 
is the only ferromagnet within this series. Especially 
CeNiSn, CePdSn, and CePtSn have been thoroughly 
investigated.* All these ternary stannides exhibit only 
one crystallographic cerium site in their structures. In 
that view, CeNiSn, CeRhSn, and CelrSn can be consid- 
ered as homogeneous intermediate- valent systems. This 
is different in the recently reported stannide CeRuSn^^ 
which is the first CeTSn compound, which adopts a su- 
perstructure (at room temperature) with two crystallo- 
graphically independent cerium sites, one trivalent and 
one intermediate-valent one. Below room temperature 
the structure becomes modulated and these modulations 
have clear consequences on the temperature dependence 
of the magnetic susceptibility and the specific heat:i£ 
Further, the system is found to order antiferromagneti- 
cally at = 3.0(2) K. Due to the peculiar structural be- 
havior and the course of the physical properties we were 
interested in the electronic structure and chemical bond- 
ing properties of CeRuSn in comparison to intermediate- 
valent CeRhSn. Our investigation is carried out in the 
framework of density functional theoretical (DFT) -'^^d^ 
and uses the scalar-relativistic implementation of the 
augmented spherical wave (ASW) methodj ^'^d^d^ 



II. CRYSTAL STRUCTURES 

The structures of CeRhSn^ and CeRuSn^ are presented 
in Fig. [1] together with the coordination polyhedra of 
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FIG. 1: The crystal structures of CeRhSn (space group 
P62m) and CeRuSn (space group C2/m) as projections along 
the short unit cell axes (bottom). The cerium, rhodium 
(ruthenium), and tin atoms are drawn as medium grey, black, 
and open circles, respectively. The trigonal prismatic units 
in CeRhSn and the three-dimensional [RuSn] network of 
CeRuSn are emphasized. The coordination polyhedra of the 
different cerium sites are presented at the top of the drawing 
together with relevant interatomic distances in units of pm. 

the cerium atoms. CeRhSn (hexagonal ZrNiAl-type) 
has only one crystallographically independent cerium 
site with five rhodium, six tin, and six cerium atoms 
in the coordination shell. The cerium coordinations of 
the two crystallographically independent cerium sites in 
CeRuSn (new monoclinic type, space group C2/m) have 
the same topology; however, small distortions lead to 
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drastically different interatomic distances. Five ruthe- 
nium, seven tin, and four cerium atoms are within the 
coordination sphere of both cerium atoms. The striking 
structural features in both ternary stannides are short 
cerium-transition metal distances, which are 304-309 pm 
Ce-Rh in CeRhSn, 233-246 pm Cel-Ru and 288-291 pm 
Ce2-Ru in CeRuSn, close to the sums of the covalent 
radii (290pm Ce + Rh and 289pm Ce -I- Ru),^^ indi- 
cating strong Ce-Rh and Ce-Ru bonding. The shortest 
interatomic distances in both structures occur for Rh-Sn 
(277-285 pm) and Ru-Sn (265-290 pm), which compare 
well with the sums of the covalent radii (265 pm Rh -f- 
Sn and 264pm Ru + Sn). Thus, the Ce-T and T-Sn 
interactions play the dominant role in chemical bond- 
ing in CeRhSn and CeRuSn. In the present work, we 
perform a comparative study of the chemical bonding 
in these two compounds as based on density functional 
ab initio calculations. In doing so, we take into account 
previous calculations as well as experimental XPS stud- 
ies on CeRhSni2iiilii2iii^i2£ Our calculations are based on 
the single-crystal data given in Ref^i^. 



III. THEORETICAL FRAMEWORK 

As a matter of fact, the degree of delocalization of 
the cerium 4/ states depends on the applied pressure as 
well as on the environment in the crystal. In electronic 
structure calculations this delicate situation is addressed 
through various approaches treating the 4/ states either 
as atomic like core states or as part of the valence basis 
set. This duality was experimentally evidenced in a com- 
bined analysis of /iSR (muon spin relaxation) and neu- 
tron experiments on cerium intermetallic systems, which 
reveals the existence of magnetic excitations due to both 
conduction electrons at the Fermi level and well local- 
ized /-electrons)^ Regarding the local environment, the 
quantum mixing (hybridization) of the 4/ states with 
those of the ligand states can have large effects as well. 
This involves chemical bonding properties, which depend 
on the crystal structure as it is illustrated by the com- 
pounds under study. 



A. Computational method 

For the electronic structure calculations we used the 
augmented spherical wave (ASW) method in its scalar- 
relativistic implementationi^iiiii^. In the ASW method 
the wave function is expanded in atom-centered aug- 
mented spherical waves, which are Hankel functions and 
numerical solutions of Schrodinger's equation, respec- 
tively, outside and inside the so-called augmentation 
spheres. In order to optimize the basis set, additional 
augmented spherical waves were placed at carefully se- 
lected interstitial sites. The choice of these sites as well 
as the augmentation radii were automatically determined 
using the sphere-geometry optimization algorithm. All 



valence states, including the Ce 4/ orbitals, were treated 
as band states. In the minimal ASW basis set, we chose 
the outermost shells to represent the valence states us- 
ing partial waves up to / = 4 for Ce as well as Z = 3 
for Ru, Rh and Sn. The completeness of the valence 
basis set was checked for charge convergence. The self- 
consistent field calculations are run to a convergence of 
AQ — lO^^Ryd'^'^ and the accuracy of the method is in 
the range of about 1 meV regarding energy differences. 



B. Spin-dependent calculations 

Due to the intermediate valent and trivalent natures 
identified for cerium within CeRuSn we carried out spin- 
degenerate non-magnetic as well as spin-polarized calcu- 
lations. Assuming a non-magnetic configuration means 
that the spin degeneracy is enforced for all species. Of 
course, such a configuration should not be confused with 
a paramagnet, which could be simulated either by a su- 
percell calculation with random spin orientations or by 
calling for disordered local moment approaches based 
on the coherent potential CPA approximation^^ or the 
LDA+DMFT scheme-^^^ 

Subsequent spin-polarized calculations with different 
initial spin populations can lead at self-consistency ci- 
ther to finite or zero local moments within an implicit 
long range ferromagnetic order. They allow to confirm 
trends established within the mean-field analysis. Anti- 
ferromagnetic (AF) calculations can be carried out to test 
for the AF ground state. This can be done for instance 
by enforcing an up-spin orientation for half of the crystal- 
lographic sites and a down-spin orientation for the other 
half thus creating two magnetic substructures. Such a 
procedure is followed here to identify the AF ground state 
of CeRuSn, which contains four formula units per cell. 



C. Assessment of chemical bonding properties 

To extract more information about the nature of the 
interactions between the atomic constituents from elec- 
tronic structure calculations, the crystal orbital over- 
lap population (COOP)^^ or the crystal orbital Hamil- 
tonian population (COHP)— may be employed. The 
latter has been already used to study chemical bond- 
ing in CeRhSn.^ While both the COOP and COHP ap- 
proaches provide a qualitative description of the bond- 
ing, nonbonding, and antibonding interactions between 
two atoms, the COOP description in some cases exag- 
gerates the magnitude of antibonding states. A slight 
refinement was recently proposed in form of the so-called 
covalent bond energy Ecov, which combines the COHP 
and COOP to calculate quantities independent of the 
particular choice of the potential zero.— In the present 
work this covalent bond energy was used for the chemi- 
cal bonding analysis. In the plots, negative, positive and 
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zero magnitudes of Ecov are indicative of bonding, anti- 
bonding, and nonbonding interactions respectively. 



IV. RESULTS FOR NON MAGNETIC 
CONFIGURATIONS 

A. Site-projected density of states 

The site-projected DOS (PDOS) for the two ternary 
stannides CeRhSn and CcRuSn arc given in Fig. [2l In all 
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FIG. 2: Non magnetic site projected DOS of CeRhSn (a) and 
CeRuSn (b) 



panels the Fermi level [Ep) is taken as zero energy. The 
cerium DOS are seen to prevail through the large peak 
around Ep mainly due to the 4/ states. The PDOS of 
CeRhSn show similarities with the plots obtained from 
the previous calclations^'^^, mainly for the itinerant part 
in the energy range from -4eV to Ep, for which Ce XPS 
spectra were measured."!^ However, major differences oc- 
cur for CeRuSn, which contains two cerium sites. Broad- 
ened Cel PDOS are found mainly above the Fermi level 
with a small contribution at Ep. In contrast, the Ce2 



PDOS are localized at Ep, which crosses the lower en- 
ergy part of the Ce 4/ states, and behave similarly to 
the Ce 4/ states in CeRhSn. As compared to the lat- 
ter, the larger splitting of the PDOS, mainly Cel as well 
as Ce2 is due to the lower symmetry of the monoclinic 
structure of the Ru-based system. There is a non negligi- 
ble contribution from Ce itinerant states below Ep which 
ensure for the chemical bonding through the hybridiza- 
tion with the transition metal (Ru, Rh) Ad and Sn 5s 
and 5p states. Due to the large filling of their d-states, 
the Ru and Rh PDOS are found below Ep, completely 
within the valence band (VB). The different PDOS show 
similar shapes within the VB, mainly in the range from 
-5eV up to Ep, which is indicative of the hybridization 
of the valence states involving itinerant Ce states below 
Ep. We will return to this issue while discussing the 
chemical bonding in terms of the covalent bond energy 
Ecov below. 



B. Analysis of the DOS within Stoner theory 

In as far as the Ce 4/ states were treated as band states 
in the framework of our calculations the Stoner theory 
of band ferromagnetism^ can be applied to address the 
spin polarization at the different cerium sites. At zero 
temperature (ground state) one can express the total en- 
ergy of the spin system resulting from the exchange and 
kinetic energies counted from a non-magnetic state as 
A[l—In{Ep)\. Here / is the Stoner exchange- 



correlation integral, which is an atomic quantity derived 
from spin-polarized scalar-relativistic calculations^, and 
n{Ep) is the DOS at Ep in the non-magnetic state. 
From this expression, the product In{Ep) provides a 
criterion for the stability of the spin system: A spin- 
polarized configuration (unequal spin occupation) will 
be more favourable if In(Ep) > 1. The system then 
stabilizes through a gain of exchange energy. From the 
calculations we have for n{Ep) for Cel, Ce2 in CeRuSn 
and Ce in the rhodium based stannide 11, 98, and 39 
states per Ryd, respectively. To analyze these results we 
use the Stoner integral for Ce as obtained from former 
scalar-relativistic calculations^^: /(Ce4/) ~ 0.02 Ryd. 
The resulting Stoner products In{Ep) for Cel, Ce2 and 
Ce are, respectively, 0.22, 1.96 and 0.78. The Stoner 
criterion is obeyed only for Ce2, which should carry a 
magnetic moment within CeRuSn when spin polarized 
calculations are carried out. This agrees with its triva- 
lent character versus the intermediate valent character of 
Cel. Still, the rather large magnitude of In{Ep) for Ce 
in CcRhSn points to a tendency towards a magnetic in- 
stability. However, according to the experiment CeRhSn 
is an intermediate- valence ternary stannide. For this rea- 
son, no ordered magnetic moment should develop when 
spin polarization is allowed for. This is further checked 
with the spin-polarized calculations outlined in the fol- 
lowing section. 
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C. Covalent bond energy Ecov 



2. CeRuSn 



Chemical bonding properties can be already addressed 
on the basis of the spin-degenerate calculations. This 
is due to the fact that the spin-polarized bands, to a 
large degree, result from the spin-degenerate bands by 
a rigid energy shift. Negative, positive and zero values 
of Ecov indicate bonding, antibonding and nonbonding 
states, respectively. 



1. CeRhSn 

The covalent bond energies Ecov for the different pairs 
of orbitals within CeRhSn are shown in Fig. [31 Here, we 
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FIG. 3: Chemical bonding: non magnetic Ecov for CeRhSn. 



have combined contributions of both Rh sites in order to 
keep the representation simple. Again, energies are re- 
ferred to the Fermi level. Note that the Ecov intensities 
along the ordinate are unitless and should be considered 
at a qualitative level. From the bottom of the valence 
band (VB) up to the Fermi level the dominant inter- 
actions result from the Ce-Rh and Rh-Sn bonds, which 
are found strongly bonding up to -2 eV. While the Ce-Rh 
bond has a rather nonbonding character up to Ep, it only 
starts to be antibonding above Ep within the empty con- 
duction band (CB). One can also notice an antibonding 
peak above Ep, which can also be observed for Ce2-Ru2 
in CeRuSn (see Fig. [4]) as it will be shown below. To the 
contrary, antibonding Rh-Sn are observed early within 
the VB, which destablizes the system. To conclude, the 
Ce-Rh bond is evidently the stabilizing contribution, fol- 
lowed by the Ce-Sn one, which is of smaller magnitude. 
It is interesting to note that, although the definition of 
the covalent bond energy Ecov is somewhat different from 
that of the COHP, the same bonding trends of CeRhSn 
are obtained from the calculations of Schmidt et al. using 
the latter 1^ 



Due to the superstructure of this system it is useful to 
study the interactions within two subsystems separately. 
The result is presented in Fig. 01 There, only one atom 
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FIG. 4: Chemical bonding: non magnetic Ecov for CeRuSn 
in two sub cells. 



of each kind is involved. As a consequence, the abso- 
lute Ecov values are smaller than in Fig. [3l The major 
interactions in the two subsystems occur clearly for Ce- 
Ru, which is bonding throughout the VB but with twice 
larger intensity for Cel than for Ce2. This follows from 
the relative distances whereby 233 < dcei-R.u < 246 pm 
while 288 < dcc2-Ru < 291pm. This is concomitant 
with the different behavior of the two cerium sites as 
discussed above. The larger Ce2-Ru separation leads to 
a stronger localization of Ce2 states, which fact is fa- 
vorable for the onset of an atomic magnetic moment. 
We also note that the Ce-Rh separation is even larger 
{dce-Rh ~ 307 pm) in CeRhSn. The antibonding peak 
above the Fermi level just like in the Rh based system 
is in line with the magnetic instability of Ce2. Finally, 
albeit with a much smaller magnitude, the Ce-Sn interac- 
tion is bonding throughout the VB and should contribute 



5 



to the bonding within the system. 

V. RESULTS OF THE SPIN POLARIZED 
CONFIGURATIONS 

Spin-polarized calculations for the magnetic structures 
were carried out by initially allowing for two different 
spin occupations, then self-consistently converging the 
charges and the magnetic moments. We assume firstly a 
hypothetic ferromagnetic configuration without any con- 
straint on the spins. Then antiferromagnetic computa- 
tions were carried out for CeRuSn in order to identify 
the ground state from energy differences. 
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A. CeRhSn 

Independent of the experimental finding of the inter- 
mediate valent character of Ce within CeRhSn, it was 
suggested from our spin-degenerate calculations as well 
as from Ref.J-i that the system could be on the verge of 
a magnetic instability. However, we could not identify a 
finite magnetic moment on Ce from self-consistent com- 
putations at high precision sampling of k-space in full 
agreement with the experiments. 

B. CeRuSn 

Contrary to the former system, a stable magnetic 
solution was identified at high precision BZ sampling. 
The energy difference favoring the ferromagnetic state 
is AE — 7.14 meV per formula unit. The spin-only 
moments were M(Cel) = O.OOS^b, M(Ce2) = 0.4Vs, 
Af (Rul) = -0.015^3 and M(Ru2) = -OMl^lB■ While 
the vanishingly small moment of Cel resembles the re- 
sults of CeRhSn, the finite moment on Ce2 is in agree- 
ment with the experimental finding as to the trivalent 
behavior of cerium. The magnetic moments on Rul and 
Ru2 are of induced nature through the hybridization 
of their states with the respective Cel and Ce2 states. 
These results are illustrated at Fig. [5] showing the site 
and spin projected density of states of CeRuSn. Low ly- 
ing Sn 5s states are not shown and the energy window 
from —6 to 6 eV shows the different behavior of the differ- 
ent Ce sites, whereby the exchange splitting is observed 
only for Ce2. The main bonding characteristics follow 
the disucssion above of the non magnetic DOS. 

Spin-orbit coupling effects can be large in Ce based 
magnetic systems. This is because the localized character 
of the 4/ wave function leads to the formation of orbital 
moments. We use the orbital field (OR) scheme intro- 
duced by Brooks^ as well as Sandratskii and Kiibler— , 
which helped to account for the experimental moments 
within formerly studied Ce intermetallic system o^'^'^'* . 
The trend is that the magnitude of the orbital moment 
(L) is close to that of the spin-only (SO) moment but with 



FIG. 5: Site and spin projected DOS of CeRuSn in the ferro- 
magnetic hypothetic state. 



opposite signs, in agreement with Hund's 3rd rule. The 
L moment of cerium stems from a 4/(Ce2) occupation of 
~ 1.3 electrons whose orbital moment 2.5/iB) comes 
close to the one of an atomic orbital, namely 3fj,B as ex- 
pected from Hund's 2nd rule. This reflects an atomic-like 
character of the 4/ (Ce) subshell. A resulting ordered LS 
moment of cerium of ~ 2/LtB is then obtained; it should be 
confronted with experimental magnitudes of magnetiza- 
tion from neutron diffraction when the results are made 
available. 

Lastly, in order to check for the ground state, AF cal- 
culations were carried out. At self consistency, AE — 
— 2.3mcV in favor of antiferromagnetic ordering was ob- 
tained thus pointing to an antiferromagnetic ground state 
in agreement with the experiment. The moment carried 
by Ce2 is ±0.39/iB, slightly smaller than in the ferromag- 
netic configuration. This is likely due to the decrease of 
symmetry when the two magnetic substructures were ac- 
counted for within the unit cell. 



VI. CONCLUSION 

In this work we have undertaken theoretical investi- 
gations of the two ternary equiatomic cerium stannides 
CeRhSn and CeRuSn as based on accurate X-ray deter- 
minations. Calculated densities of states and chemical 
bonding properties for CeRhSn agree very well with for- 
mer theoretical results. CeRuSn is strongly infiuenced 
by the peculiar feature of a complex interplay of both in- 
termediate valent and trivalent Ce atoms. The exciting 
properties of this compound were addressed both by spin- 
degenerate and spin-polarized density functional theory 
based calculations. Analysis of the electronic structures 
and of the chemical bonding reveal different types of 
chemical bonds due to the nature of the Ce site and its 
environment with Ru, Sn as well as Ce ligands. As a 
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consequence, a local magnetic moment is expected only 
for the cerium site with trivalent character. Indeed, spin- 
polarized calculations lead to a nearly vanishing magnetic 
moment at the intermediate-valent Cel site and a finite 
moment shows up only at the trivalent Ce2 site. In agree- 
ment with experimental data, the Ce2 moments give rise 
to an antiferromagnetic ground state. 
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